What is claimed is: 



1 . A method of obtaining magnetic resonance signals with signal separation for at least two 
chemical species in a heterogeneous magnetic field, comprising the steps of: 

a) obtaining first magnetic resonance signals from pixels in an object having at least two 
chemical species using a first repetition time and a first echo time, 

b) obtaining at least second and third magnetic resonant signals from the pixels using 

second and third echo times, 

c) determining a signal estimate for each pixel by combining all signals for the pixel 
using a linear least squares fitting to decompose the chemical species, assuming a first value of 

field heterogeneity (M, 

d) calculating a first error to the field heterogeneity, 

e) repeating step c) using the first value of field heterogeneity and the error from step d), 

f) repeating step d) to calculate a second error to the field heterogeneity, and 

g) updating the value of field heterogeneity and repeating steps c) and d) until an 
acceptable error is calculated. 

2. The method as defined by claim lwherein M chemical species are present and step b) 
includes obtaining at least M+l magnetic resonance signals for each pixel. 

3 . The method as defined by claim 1 wherein fat and water are two chemical species and step 
b) includes obtaining three magnetic resonance signals. 

4. The method as defined by claim 1 wherein in step b) the magnetic resonance signals at 
time, n=l to N, for species, j =1 to M, having real, R, and imaginary, I, parts is given by: 

and a least squares fitting of all signals is given by: 

p = (A r A)-'A r S 
where A is a known matrix for M species. 
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5. The method as defined by claim 4 wherein: 
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where c? = cosfca^O' d ? = ™fa*f/J*)> s" =2^/,(- p' w - p R f d„- p' f c m ) and 
5 g'„ = 2;rr„ (p* + p;c n - ^„ ) are the matrix elements. 

6. The method as defined by claim 4 wherein error to the field heterogeneity is given by: 

y = (B r B)-'B r i 
where B is a matrix given by: 
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where 

s = \a ¥ Ap* Ap 1 Ap R 2 Ap[ ... Ap K M Ap 1 M \,gl=2nt n Y\-p*d jn -p)c jn )K 

1 1 7=1 

g'j,=2xt n ±(p5c Jn -p' j d Jn l 

15 7. The method as defined by claim 6 wherein step a) includes obtaining signals from a 
single coil. 



STFUP152/S03-130 



Page 12 



8. The method as defined by claim 6 wherein step a) includes obtaining signals from a 
plurality of coils and steps b) through g) are performed for signals from each coil, and further 

including the step of: 

h) combining field heterogeneity as determined from signals for each coil. 

9. The method as defined by claim 8 wherein field heterogeneity is determined by weighting 
contributions from each coil. 

1 0. The method as defined by claim 9 wherein the weighting contribution from each coil is a 
function of the square of the magnitude of the image contributed by that coil. 

11. The method as defined by claim 1 0 where for each pixel, at position r, the combined field 
map is: 




where P coils collect P independent images. 

12. The method as defined by claim 1 1 wherein the combined field heterogeneity from step 
h) is smoothed by passing through a low pass filter. 

13. The method as defined by claim 1 1 wherein M images of each species are obtained using 
signals from each of M coils and the combined field heterogeneity, and then combining the M 
images using a square root of the sum of the squares of the images. 

14. The method as defined by claim 1 wherein step a) includes obtaining signals from a 
plurality of coils and steps b) through g) are performed for signals from each coil, and further 

including the step of: 

h) combining field heterogeneity as determined from signals for each coil. 

15. The method as defined by claim 14 wherein field heterogeneity is determined by 
weighting contributions from each coil. 
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16. The method as defined by claim 15 wherein the weighting contribution from each coil is a 
function of the square of the magnitude of the image contributed by that coil. 

17. The method as defined by claim 16 where for each pixel, at position r, the combined field 
5 map is: 




where P coils collect P independent images. 

18. The method as defined by claim 17 wherein the combined field heterogeneity from step h) is 
1 0 smoothed by passing through a low pass filter. 

19. The method as defined by claim 17 wherein M images of each species is obtained using 
signals from each of M coils and the combined field heterogeneity, and then combining the M 
images using a square root of the sum of the squares of the images. 

15 
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